Optimal population-specific HLA imputation with dimension reduction.

Human genomics has quickly evolved, powering genome-wide association studies (GWASs). SNP-based GWASs cannot capture the intense polymorphism of HLA genes, highly associated with disease susceptibility. There are methods to statistically impute HLA genotypes from SNP-genotypes data, but lack of diversity in reference panels hinders their performance. We evaluated the accuracy of the 1000 Genomes data as a reference panel for imputing HLA from admixed individuals of African and European ancestries, focusing on (a) the full dataset, (b) 10 replications from 6 populations, and (c) 19 conditions for the custom reference panels. The full dataset outperformed smaller models, with a good F1-score of 0.66 for HLA-B. However, custom models outperformed the multiethnic or population models of similar size (F1-scores up to 0.53, against up to 0.42). We demonstrated the importance of using genetically specific models for imputing populations, which are currently underrepresented in public datasets, opening the door to HLA imputation for every genetic population.

18th International HLA and Immunogenetics Workshop: Report on the SNP-HLA Reference Consortium (SHLARC) component.

The SNP-HLA Reference Consortium (SHLARC), a component of the 18th International HLA and Immunogenetics Workshop, is aimed at collecting diverse and extensive human leukocyte antigen (HLA) data to create custom reference panels and enhance HLA imputation techniques. Genome-wide association studies (GWAS) have significantly contributed to identifying genetic associations with various diseases. The HLA genomic region has emerged as the top locus in GWAS, particularly in immune-related disorders. However, the limited information provided by single nucleotide polymorphisms (SNPs), the hallmark of GWAS, poses challenges, especially in the HLA region, where strong linkage disequilibrium (LD) spans several megabases. HLA imputation techniques have been developed using statistical inference in response to these challenges. These techniques enable the prediction of HLA alleles from genotyped GWAS SNPs. Here we present the SHLARC activities, a collaborative effort to create extensive, and multi-ethnic reference panels to enhance HLA imputation accuracy.

Mitochondrial DNA haplogroups and variants predispose to chagas disease cardiomyopathy

Cardiomyopathies are major causes of heart failure. Chagas disease (CD) is caused by the parasite Trypanosoma cruzi, and it is endemic in Central and South America. Thirty percent of cases evolve into chronic chagas cardiomyopathy (CCC), which has worse prognosis as compared with other cardiomyopathies. In vivo bioenergetic analysis and ex vivo proteomic analysis of myocardial tissues highlighted worse mitochondrial dysfunction in CCC, and previous studies identified nuclear-encoded mitochondrial gene variants segregating with CCC. Here, we assessed the role of the mitochondrial genome through mtDNA copy number variations and mtDNA haplotyping and sequencing from heart or blood tissues of severe, moderate CCC and asymptomatic/indeterminate Chagas disease as well as healthy controls as an attempt to help decipher mitochondrial-intrinsic genetic involvement in Chagas disease development. We have found that the mtDNA copy number was significantly lower in CCC than in heart tissue from healthy individuals, while blood mtDNA content was similar among asymptomatic Chagas disease, moderate, and severe CCC patients. An MtDNA haplogrouping study has indicated that African haplogroups were over represented in the Chagas subject groups in comparison with healthy Brazilian individuals. The European lineage is associated with protection against cardiomyopathy and the macro haplogroup H is associated with increased risk towards CCC. Using mitochondria DNA sequencing, 84 mtDNA-encoded protein sequence pathogenic variants were associated with CCC. Among them, two variants were associated to left ventricular non-compaction and two to hypertrophic cardiomyopathy. The finding that mitochondrial protein-coding SNPs and mitochondrial haplogroups associate with risk of evolving to CCC is consistent with a key role of mitochondrial DNA in the development of chronic chagas disease cardiomyopathy.

Association between APOE-ε4 allele and cognitive function is mediated by Alzheimer's disease pathology: a population-based autopsy study in an admixed sample.

BACKGROUND: Apolipoprotein E ε4 allele (APOE-ε4) is the main genetic risk factor for late-onset Alzheimer's disease (AD) and may impact cognitive function also via other neuropathological lesions. However, there is limited evidence available from diverse populations, as APOE associations with dementia seem to differ by race. Therefore, we aimed to evaluate the pathways linking APOE-ε4 to cognitive abilities through AD and non-AD neuropathology in an autopsy study with an admixed sample. METHODS: Neuropathological lesions were evaluated following international criteria using immunohistochemistry. Participants were classified into APOE-ε4 carriers (at least one ε4 allele) and non-carriers. Cognitive abilities were evaluated by the Clinical Dementia Rating Scale sum of boxes. Mediation analyses were conducted to assess the indirect association of APOE-ε4 with cognition through AD-pathology, lacunar infarcts, hyaline arteriosclerosis, cerebral amyloid angiopathy (CAA), Lewy body disease (LBD), and TAR DNA-binding protein 43 (TDP-43). RESULTS: We included 648 participants (mean age 75 ± 12 years old, mean education 4.4 ± 3.7 years, 52% women, 69% White, and 28% APOE-ε4 carriers). The association between APOE-ε4 and cognitive abilities was mediated by neurofibrillary tangles (ß = 0.88, 95% CI = 0.45; 1.38, p < 0.001) and neuritic plaques (ß = 1.36, 95% CI = 0.86; 1.96, p < 0.001). Lacunar infarcts, hyaline arteriosclerosis, CAA, LBD, and TDP-43 were not mediators in the pathway from APOE-ε4 to cognition. CONCLUSION: The association between APOE-ε4 and cognitive abilities was partially mediated by AD-pathology. On the other hand, cerebrovascular lesions and other neurodegenerative diseases did not mediate the association between APOE-ε4 and cognition.

Apolipoprotein E ε2 allele is associated with lower risk of carotid artery obstruction in a population-based autopsy study.

INTRODUCTION: Apolipoprotein E (APOE) ε4 allele has been associated with higher carotid atherosclerosis risk, while the APOE-ε2 seems to decrease this risk. Data from autopsy studies, where carotid arteries can be evaluated in their full extension, is scarce. Therefore, we investigated the association between APOE alleles and direct morphometric measurements of carotid atherosclerosis in an autopsy study with an admixed sample. METHODS: We measured the intima-media thickness (IMT) and stenosis of the common (CCA) and internal carotid (ICA) arteries. The APOE polymorphisms were determined by real-time polymerase chain reaction. Participants were classified into three groups according to the APOE alleles (ε2, ε3, and ε4). We evaluated the association between APOE groups and carotid atherosclerosis using adjusted regression models and included interaction terms of APOE alleles with age, sex, and race. RESULTS: We evaluated 1,850 carotid artery samples from 185 participants (mean age=75±12 years old, 55% female, and 71% White). The APOE-ε2 group (n=17) had a lower carotid obstruction and a lower number of severe stenoses (≥ 70%). Having at least one ε4 allele (n=51) was not associated with carotid atherosclerosis. APOE alleles were also not associated with carotid IMT. Age, sex, and race did not modify these relationships. CONCLUSION: APOE-ε2 carriers had a lower percentage of carotid obstruction and less severe stenosis. APOE-ε4 was not related to a higher risk of carotid atherosclerosis in this cross-sectional population-based autopsy study.

Identification of pathogenic variants in the brazilian cohort with familial hypercholesterolemia using exon-targeted gene sequencing

Familial hypercholesterolemia (FH) is a monogenic disease characterized by high plasma low-density lipoprotein cholesterol (LDL-c) levels and increased risk of premature atherosclerotic cardiovascular disease. Mutations in FH-related genes account for 40% of FH cases worldwide. In this study, we aimed to assess the pathogenic variants in FH-related genes in the Brazilian FH cohort FHBGEP using exon-targeted gene sequencing (ETGS) strategy. FH patients (n = 210) were enrolled at five clinical sites and peripheral blood samples were obtained for laboratory testing and genomic DNA extraction. ETGS was performed using MiSeq platform (Illumina). To identify deleterious variants in LDLR, APOB, PCSK9, and LDLRAP1, the long-reads were subjected to Burrows-Wheeler Aligner (BWA) for alignment and mapping, followed by variant calling using Genome Analysis Toolkit (GATK) and ANNOVAR for variant annotation. The variants were further filtered using in-house custom scripts and classified according to the American College Medical Genetics and Genomics (ACMG) guidelines. A total of 174 variants were identified including 85 missense, 3 stop-gain, 9 splice-site, 6 InDel, and 71 in regulatory regions (3'UTR and 5'UTR). Fifty-two patients (24.7%) had 30 known pathogenic or likely pathogenic variants in FH-related genes according to the American College Medical and Genetics and Genomics guidelines. Fifty-three known variants were classified as benign, or likely benign and 87 known variants have shown uncertain significance. Four novel variants were discovered and classified as such due to their absence in existing databases. In conclusion, ETGS and in silico prediction studies are useful tools for screening deleterious variants and identification of novel variants in FH-related genes, they also contribute to the molecular diagnosis in the FHBGEP cohort.

Characterization of pharmacogenomic variants in a Brazilian admixed cohort of elderly individuals based on whole-genome sequencing data.

Introduction: Research in the field of pharmacogenomics (PGx) aims to identify genetic variants that modulate response to drugs, through alterations in their pharmacokinetics (PK) or pharmacodynamics (PD). The distribution of PGx variants differs considerably among populations, and whole-genome sequencing (WGS) plays a major role as a comprehensive approach to detect both common and rare variants. This study evaluated the frequency of PGx markers in the context of the Brazilian population, using data from a population-based admixed cohort from Sao Paulo, Brazil, which includes variants from WGS of 1,171 unrelated, elderly individuals. Methods: The Stargazer tool was used to call star alleles and structural variants (SVs) from 38 pharmacogenes. Clinically relevant variants were investigated, and the predicted drug response phenotype was analyzed in combination with the medication record to assess individuals potentially at high-risk of gene-drug interaction. Results: In total, 352 unique star alleles or haplotypes were observed, of which 255 and 199 had a frequency < 0.05 and < 0.01, respectively. For star alleles with frequency > 5% (n = 97), decreased, loss-of-function and unknown function accounted for 13.4%, 8.2% and 27.8% of alleles or haplotypes, respectively. Structural variants (SVs) were identified in 35 genes for at least one individual, and occurred with frequencies >5% for CYP2D6, CYP2A6, GSTM1, and UGT2B17. Overall 98.0% of the individuals carried at least one high risk genotype-predicted phenotype in pharmacogenes with PharmGKB level of evidence 1A for drug interaction. The Electronic Health Record (EHR) Priority Result Notation and the cohort medication registry were combined to assess high-risk gene-drug interactions. In general, 42.0% of the cohort used at least one PharmGKB evidence level 1A drug, and 18.9% of individuals who used PharmGKB evidence level 1A drugs had a genotype-predicted phenotype of high-risk gene-drug interaction. Conclusion: This study described the applicability of next-generation sequencing (NGS) techniques for translating PGx variants into clinically relevant phenotypes on a large scale in the Brazilian population and explores the feasibility of systematic adoption of PGx testing in Brazil.

Corrigendum: Follow-up of young adult monozygotic twins after simultaneous critical coronavirus disease 2019: a case report.

[This corrects the article DOI: 10.3389/fmed.2022.1008585.].

Identification of pathogenic variants in the Brazilian cohort with Familial hypercholesterolemia using exon-targeted gene sequencing.

Familial hypercholesterolemia (FH) is a monogenic disease characterized by high plasma low-density lipoprotein cholesterol (LDL-c) levels and increased risk of premature atherosclerotic cardiovascular disease. Mutations in FH-related genes account for 40% of FH cases worldwide. In this study, we aimed to assess the pathogenic variants in FH-related genes in the Brazilian FH cohort FHBGEP using exon-targeted gene sequencing (ETGS) strategy. FH patients (n = 210) were enrolled at five clinical sites and peripheral blood samples were obtained for laboratory testing and genomic DNA extraction. ETGS was performed using MiSeq platform (Illumina). To identify deleterious variants in LDLR, APOB, PCSK9, and LDLRAP1, the long-reads were subjected to Burrows-Wheeler Aligner (BWA) for alignment and mapping, followed by variant calling using Genome Analysis Toolkit (GATK) and ANNOVAR for variant annotation. The variants were further filtered using in-house custom scripts and classified according to the American College Medical Genetics and Genomics (ACMG) guidelines. A total of 174 variants were identified including 85 missense, 3 stop-gain, 9 splice-site, 6 InDel, and 71 in regulatory regions (3'UTR and 5'UTR). Fifty-two patients (24.7%) had 30 known pathogenic or likely pathogenic variants in FH-related genes according to the American College Medical and Genetics and Genomics guidelines. Fifty-three known variants were classified as benign, or likely benign and 87 known variants have shown uncertain significance. Four novel variants were discovered and classified as such due to their absence in existing databases. In conclusion, ETGS and in silico prediction studies are useful tools for screening deleterious variants and identification of novel variants in FH-related genes, they also contribute to the molecular diagnosis in the FHBGEP cohort.

Dynamics of the personalities of PCSK9 on missense variants (rs505151 and rs562556) from elderly cohort studies in Brazil.

The Proprotein Convertase Subtilisin/Kexin Type 9 (PCSK9) promotes the degradation of the low-density lipoprotein receptors (LDLR). Gain-of-function (GOF) variants of PCSK9 significantly affects lipid metabolism leading to coronary artery disease (CAD), owing to the raising the plasma low-density lipoprotein (LDL). Considering the public health matter, large-scale genomic studies have been conducted worldwide to provide the genetic architecture of populations for the implementation of precision medicine actions. Nevertheless, despite the advances in genomic studies, non-European populations are still underrepresented in public genomic data banks. Despite this, we found two high-frequency variants (rs505151 and rs562556) in the ABraOM databank (Brazilian genomic variants) from a cohort SABE study conducted in the largest city of Brazil, São Paulo. Here, we assessed the structural and dynamical features of these variants against WT through a molecular dynamics study. We sought fundamental dynamical interdomain relations through Perturb Response Scanning (PRS) and we found an interesting change of dynamical relation between prodomain and Cysteine-Histidine-Rich-Domain (CHRD) in the variants. The results highlight the pivotal role of prodomain in the PCSK9 dynamic and the implications for the development of new drugs depending on patient group genotype.

Immunogenetics of HLA-B: SNP, allele, and haplotype diversity in populations from different continents and ancestry backgrounds.

HLA-B is among the most variable gene in the human genome. This gene encodes a key molecule for antigen presentation to CD8+ T lymphocytes and NK cell modulation. Despite the myriad of studies evaluating its coding region (with an emphasis on exons 2 and 3), few studies evaluated introns and regulatory sequences in real population samples. Thus, HLA-B variability is probably underestimated. We applied a bioinformatics pipeline tailored for HLA genes on 5347 samples from 80 different populations, which includes more than 1000 admixed Brazilians, to evaluate the HLA-B variability (SNPs, indels, MNPs, alleles, and haplotypes) in exons, introns, and regulatory regions. We observed 610 variable sites throughout HLA-B; the most frequent variants are shared worldwide. However, the haplotype distribution is geographically structured. We detected 920 full-length haplotypes (exons, introns, and untranslated regions) encoding 239 different protein sequences. HLA-B gene diversity is higher in admixed populations and Europeans while lower in African ancestry individuals. Each HLA-B allele group is associated with specific promoter sequences. This HLA-B variation resource may improve HLA imputation accuracy and disease-association studies and provide evolutionary insights regarding HLA-B genetic diversity in human populations.

Immunological evaluation of young unvaccinated patients with Turner syndrome after COVID-19.

OBJECTIVES: The X-chromosome contains the largest number of immune-related genes, which play a major role in COVID-19 symptomatology and susceptibility. Here, we had a unique opportunity to investigate, for the first time, COVID-19 outcomes in six unvaccinated young Brazilian patients with Turner syndrome (TS; 45, X0), including one case of critical illness in a child aged 10 years, to evaluate their immune response according to their genetic profile. METHODS: A serological analysis of humoral immune response against SARS-CoV-2, phenotypic characterization of antiviral responses in peripheral blood mononuclear cells after stimuli, and the production of cytotoxic cytokines of T lymphocytes and natural killer cells were performed in blood samples collected from the patients with TS during the convalescence period. Whole exome sequencing was also performed. RESULTS: Our volunteers with TS showed a delayed or insufficient humoral immune response to SARS-CoV-2 (particularly immunoglobulin G) and a decrease in interferon-γ production by cluster of differentiation (CD)4+ and CD8+ T lymphocytes after stimulation with toll-like receptors 7/8 agonists. In contrast, we observed a higher cytotoxic activity in the volunteers with TS than the volunteers without TS after phorbol myristate acetate/ionomycin stimulation, particularly granzyme B and perforin by CD8+ and natural killer cells. Interestingly, two volunteers with TS carry rare genetic variants in genes that regulate type I and III interferon immunity. CONCLUSION: Following previous reports in the literature for other conditions, our data showed that patients with TS may have an impaired immune response against SARS-CoV-2. Furthermore, other medical conditions associated with TS could make them more vulnerable to COVID-19.

Apolipoprotein E genotypes were not associated with intracranial atherosclerosis: a population-based autopsy study.

BACKGROUND: Apolipoprotein E gene (APOE) ε4 allele is associated with a higher risk of carotid atherosclerosis, but less is known about the association of APOE with intracranial atherosclerotic disease (IAD). We aimed to investigate the association of APOE alleles with IAD in a cross-sectional autopsy study. METHODS: We measured the stenosis in the 12 arteries of the Circle of Willis using postmortem morphometric measurements. The APOE polymorphism was determined by real-time polymerase chain reaction. We assessed the association between APOE polymorphism and IAD using regression models adjusted for sociodemographic and clinical variables. We also verified the modifier effect of age, sex, and race on this association. We stratified the analysis by age group to investigate the possibility of attrition bias. RESULTS: In 400 participants (mean age=73.2±12.3 years old, 51% female, and 64% White), IAD was evaluated in 4,504 artery segments. APOE-ε4 was not associated with IAD nor with the number of artery stenosis compared to non-APOE-ε4 carriers. Sociodemographic variables did not modify this relationship. Among participants older than 70 years, there was a trend towards an association between APOE allele ε4 and a lower stenosis index in the middle cerebral artery, suggesting attrition bias related to the APOE-ε4 effect on mortality. CONCLUSIONS: APOE alleles were not associated with IAD in this population-based autopsy study. Lower stenosis in older participants suggests the possibility of attrition bias.

Challenges in selecting admixture models and marker sets to infer genetic ancestry in a Brazilian admixed population.

The inference of genetic ancestry plays an increasingly prominent role in clinical, population, and forensic genetics studies. Several genotyping strategies and analytical methodologies have been developed over the last few decades to assign individuals to specific biogeographic regions. However, despite these efforts, ancestry inference in populations with a recent history of admixture, such as those in Brazil, remains a challenge. In admixed populations, proportion and components of genetic ancestry vary on different levels: (i) between populations; (ii) between individuals of the same population, and (iii) throughout the individual's genome. The present study evaluated 1171 admixed Brazilian samples to compare the genetic ancestry inferred by tri-/tetra-hybrid admixture models and evaluated different marker sets from those with small numbers of ancestry informative markers panels (AIMs), to high-density SNPs (HDSNP) and whole-genome-sequence (WGS) data. Analyses revealed greater variation in the correlation coefficient of ancestry components within and between admixed populations, especially for minority ancestral components. We also observed positive correlation between the number of markers in the AIMs panel and HDSNP/WGS. Furthermore, the greater the number of markers, the more accurate the tri-/tetra-hybrid admixture models.

The oldest unvaccinated Covid-19 survivors in South America.

BACKGROUND: Although older adults are at a high risk of severe or critical Covid-19, there are many cases of unvaccinated centenarians who had a silent infection or recovered from mild or moderate Covid-19. We studied three Brazilian supercentenarians, older than 110 years, who survived Covid-19 in 2020 before being vaccinated. RESULTS: Despite their advanced age, humoral immune response analysis showed that these individuals displayed robust levels of IgG and neutralizing antibodies (NAbs) against SARS-CoV-2. Enrichment of plasma proteins and metabolites related to innate immune response and host defense was also observed. None presented autoantibodies (auto-Abs) to type I interferon (IFN). Furthermore, these supercentenarians do not carry rare variants in genes underlying the known inborn errors of immunity, including particular inborn errors of type I IFN. CONCLUSION: These observations suggest that their Covid-19 resilience might be a combination of their genetic background and their innate and adaptive immunity.

MUC22, HLA-A, and HLA-DOB variants and COVID-19 in resilient super-agers from Brazil.

Background: Although aging correlates with a worse prognosis for Covid-19, super elderly still unvaccinated individuals presenting mild or no symptoms have been reported worldwide. Most of the reported genetic variants responsible for increased disease susceptibility are associated with immune response, involving type I IFN immunity and modulation; HLA cluster genes; inflammasome activation; genes of interleukins; and chemokines receptors. On the other hand, little is known about the resistance mechanisms against SARS-CoV-2 infection. Here, we addressed polymorphisms in the MHC region associated with Covid-19 outcome in super elderly resilient patients as compared to younger patients with a severe outcome. Methods: SARS-CoV-2 infection was confirmed by RT-PCR test. Aiming to identify candidate genes associated with host resistance, we investigated 87 individuals older than 90 years who recovered from Covid-19 with mild symptoms or who remained asymptomatic following positive test for SARS-CoV-2 as compared to 55 individuals younger than 60 years who had a severe disease or died due to Covid-19, as well as to the general elderly population from the same city. Whole-exome sequencing and an in-depth analysis of the MHC region was performed. All samples were collected in early 2020 and before the local vaccination programs started. Results: We found that the resilient super elderly group displayed a higher frequency of some missense variants in the MUC22 gene (a member of the mucins' family) as one of the strongest signals in the MHC region as compared to the severe Covid-19 group and the general elderly control population. For example, the missense variant rs62399430 at MUC22 is two times more frequent among the resilient super elderly (p = 0.00002, OR = 2.24). Conclusion: Since the pro-inflammatory basal state in the elderly may enhance the susceptibility to severe Covid-19, we hypothesized that MUC22 might play an important protective role against severe Covid-19, by reducing overactive immune responses in the senior population.

Follow-up of young adult monozygotic twins after simultaneous critical coronavirus disease 2019: A case report.

Background: The influence of the host genome on coronavirus disease 2019 (COVID-19) susceptibility and severity is supported by reports on monozygotic (MZ) twins where both were infected simultaneously with similar disease outcomes, including several who died due to the SARS-CoV-2 infection within days apart. However, successive exposures to pathogens throughout life along with other environmental factors make the immune response unique for each individual, even among MZ twins. Case presentation and methods: Here we report a case of a young adult monozygotic twin pair, who caught attention since both presented simultaneously severe COVID-19 with the need for oxygen support despite age and good health conditions. One of the twins, who spent more time hospitalized, reported symptoms of long-COVID even 7 months after infection. Immune cell profile and specific responses to SARS-CoV-2 were evaluated as well as whole exome sequencing. Conclusion: Although the MZ twin brothers shared the same genetic mutations which may be associated with their increased risk of developing severe COVID-19, their clinical progression was different, reinforcing the role of both immune response and genetics in the COVID-19 presentation and course. Besides, post-COVID syndrome was observed in one of them, corroborating an association between the duration of hospitalization and the occurrence of long-COVID symptoms.

Global and local ancestry modulate APOE association with Alzheimer's neuropathology and cognitive outcomes in an admixed sample.

Dementia is more prevalent in Blacks than in Whites, likely due to a combination of environmental and biological factors. Paradoxically, clinical studies suggest an attenuation of APOE ε4 risk of dementia in African ancestry (AFR), but a dearth of neuropathological data preclude the interpretation of the biological factors underlying these findings, including the association between APOE ε4 risk and Alzheimer's disease (AD) pathology, the most frequent cause of dementia. We investigated the interaction between African ancestry, AD-related neuropathology, APOE genotype, and functional cognition in a postmortem sample of 400 individuals with a range of AD pathology severity and lack of comorbid neuropathology from a cohort of community-dwelling, admixed Brazilians. Increasing proportions of African ancestry (AFR) correlated with a lower burden of neuritic plaques (NP). However, for individuals with a severe burden of NP and neurofibrillary tangles (NFT), AFR proportion was associated with worse Clinical Dementia Rating sum of boxes (CDR-SOB). Among APOE ε4 carriers, the association between AFR proportion and CDR-SOB disappeared. APOE local ancestry inference of a subset of 309 individuals revealed that, in APOE ε4 noncarriers, non-European APOE background correlated with lower NP burden and, also, worse cognitive outcomes than European APOE when adjusting by NP burden. Finally, APOE ε4 was associated with worse AD neuropathological burden only in a European APOE background. APOE genotype and its association with AD neuropathology and clinical pattern are highly influenced by ancestry, with AFR associated with lower NP burden and attenuated APOE ε4 risk compared to European ancestry.